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Abstract 

Repetition-based retransmission is used in Alamouti-modulation [1998] for 2x2 MIMO sys- 
tems. We propose to use instead of ordinary repetition so-called "scaled repetition" together 
with rotation. It is shown that the rotated and scaled Alamouti code has a hard-decision 
performance which is only slightly worse than that of the Golden code [2005], the best known 
2x2 space-time code. Decoding the Golden code requires an exhaustive search over all 
codewords, while our rotated and scaled Alamouti code can be decoded with an acceptable 
complexity however. 



1 Scaled-repetition Retransmission for the SISO Channel 

First we consider transmission over a single-input single-output (SISO) additive white Gaussian 
noise (AWGN) channel (see Fig. [1]), and introduce scaled-repetition retransmission. It turns out 
that scaled-repetition improves upon ordinary-repetition retransmission. 

1.1 Some information theory 
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Figure 1: The AWGN channel. 



The real- valued output y k for transmission fc = 1,2,- 

Vk = x k + n k , 



, K, see Fig. [TJ satisfies 



(1) 



where x k is the real- valued channel input for transmission k and n k is a real- valued Gaussian noise 
sample with mean E[Nk] = 0, variance E[N?] — cr 2 , which is uncorrelated with all other noise 
samples. The transmitter power is limited, i.e. we require that -E[A^] < P. It is well-known that 
an X which is Gaussian with mean and variance P achieves capacity. This basic capacity (in 
bit/transm.) equals 

C7=~log 2 (l + J). (2) 

When we retransmit (repeat) codewords , each symbol Xk from such a codeword (x%, X2, ■ ■ ■ , xk) 
is actually transmitted and received twice, i.e. Xki = Xk2 — %ki and 



Uki = Xk + n k i, and yui — Xk + nui- 



(3) 
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Figure 2: Two mappings from Xki to Xk 2 ■ On the right the scaled-repetition mapping, left the 
ordinary-repetition mapping. 



An optimal receiver can form z k = yhl Z, Vk2 = x k + ""+" fc2 . Now the variance of the noise variable 
(Nki + Nki)/'} is <t 2 /2. Therefore the repetition capacity for a single repetition in bit/transm. is 

1 2P 

a = i io g2 (i + ^). (4) 

Fig. [3J shows the basic capacity C (black line) and repetition capacity C r (blue line) as a function 
of the signal-to-noise ratio SNR which is defined as 

SNR = P/(t 2 . (5) 

It is easy to see that always C r < C. For large SNR we may write C r w C/2 + 1/4, while for small 
SNR we obtain C r m C. 



1.2 Ordinary and scaled repetition for 4-PAM 

When we use 4-PAM modulation, the channel inputs Xk assume values from Aa-pam = {— 3, — 1, +1, 
+3}, each with probability 1/4. Ordinary repetition, see ©, leads to signal points (xi,JCa) = (%, x) 
for x € ^4_pam, see the left part of Fig. [21 For this case the maximum transmission rate 
I{X;Y\,Y<2) is shown in Fig. [3] with blue asterisks. Note that this maximum transmission rate 
is slightly smaller than the corresponding capacities C r , mainly because uniform inputs are used 
instead of Gaussians. 

We can use Benelli's [3] method to improve upon ordinary-repetition retransmission, i.e. by 
modulating the retransmitted symbol differently. We could e.g. take 

x k i = x k , and x k2 = M 2 (x k ) for x k S Ai-pam, (6) 

where M 2 (a) = 2a— 5 if a > and M 2 (a) = 2a + 5 for a < 0. We call this method scaled repetition 
since we scale a symbol by a factor (2 here) and then compensate (add -5 or +5) in order to obtain a 
symbol from A.pam- This results in the signal points (x, M 2 (x)) for x € Aupam, see Fig. [2l right 
part. Also for the scaled-repetition case the maximum transmission rate I(X;Yi,Y 2 ) is shown 
in figure [3J now with red asterisks. Note that this maximum transmission rate is only slightly 
smaller than the basic capacity C. Ordinary repetition is however definitively inferior to the basic 
transmission if the SNR is not very small. 



1.3 Demodulation complexity 

Scaled repetition outperforms ordinary repetition, but also has a disadvantage. In an ordinary- 
repetition system the output = {y k i + J/fe2)/2 is simply sliced. In a system that uses scaled 
repetition we can only slice after having distinguished between two cases. More precisely note 
that Xk 2 — M 2 {xk) — 2x k — D 2 (x k ), where D 2 (a) = 5 if a > and D 2 (a) = — 5 if a < 0. Now we 
can use a sheer for y k i + 2y k2 = x k + n kl + 2(2x k - D 2 (x k ) + n k2 ) = 5x k - 2D 2 (x k ) + n k i + 2n k2 . 
Assuming that x k <E {—3,-1} we get that D 2 (x k ) = —5 and this implies that we should put a 
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Figure 3: Basic capacity C (black curve) and repetition capacity C r (blue) in bit/transm. as a 
function of SNR = P/a 2 in dB (horizontally). Also the maximum transmission rates achievable 
with 4-PAM in the ordinary-repetition case (blue *'s). In red *'s the maximum rates achievable 
using scaled-repetition mapping. 

threshold at to distinguish between —3 and — 1. Similarly assuming that Xk G { + l,+3} we 
get D-iixk) = 5 and we must slice yki + ^yk2 again with a threshold at 0. Then the best overall 
candidate Xk is found by minimizing (yki — Xk) 2 + (yk2 — AfaOcfc)) over the two candidates. 

2 Fundamental Properties for the 2x2 MIMO Channel 

2.1 Model description 
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Figure 4: Model of a 2 x 2 MIMO channel. 

Next consider a 2 x 2 MIMO channel (see Fig. [5]). Both the transmitter and the receiver use 
two antennas. The output vector (yikiDik) at transmission k relates to the corresponding input 
vector (xik,X2k) as given by 

( Vik \ _ / hn h u \ ( xik \ ( nik \ 
\ V2k ) V ft 2i ^22 ) \ %2k ) \ n 2 k ) 

where (Nik, N2k) is a pair of independent zero-mean circularly symmetric complex Gaussians, 
both having variance a 2 (per two dimensions). Noise variable pairs in different transmissions are 
independent. 
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We assume that the four channel coefficients fin, i?i2, H^x, and H22 are independent zero- 
mean circularly symmetric complex Gaussians, each having variance 1 (per two dimensions). The 
channel coefficients are chosen prior to a block of K transmissions and remain constant over that 
block. 

The complex transmitted symbols (Xki, X k 2) must satisfy a power constraint, i.e. 

E[X kl X* kl + X k2 X* k2 ] < P. (8) 



2.2 Telatar capacity 

If the channel input variables are independent zero-mean circularly symmetric complex Gaussians 
both having variance P/2, then the resulting mutual information (called Telatar capacity here, 
see [5]) isQ 

CTelatax(ff) = log 2 det(7 2 + ## f ), (9) 

where H = ( ^ n (* 12 J , i.e. the actual channel-coefficient matrix and I 2 the 2x2 identity 



matrix. Also in the 2x2 MIMO case we define the signal-to-niose ratio as 

SNR = P/ct 2 . (10) 

It can be shown (see e.g. Yao ([5J, p. 36) that for fixed R and SNR large enough Pr{C Tclatar (#) < 
R} w 7 • SNR" 4 , for some constant 7. 



2.3 Worst-case error-probabilities 

Consider M (one for each message) K x 2 code- matrices c 1; c 2 , • • • , c M resulting in a unit average 
energy code. Then Tarokh, Seshadri and Calderbank [3] showed that for large SNR 

Pr{c -> c'} « 7'(det((c' - c)(c' - c)t)- 2 SNR" 4 . (11) 

for some 7' if the rank of the difference matrices c — c' is 2, and we transmit x = \[Pc. If this 
holds for all difference matrices we say that the diversity order is 4. Therefore it makes sense to 
maximize the minimum modulus of the determinant over all code-matrix differences. 



3 Alamouti: Ordinary Repetition 

Alamouti PQ proposed a modulation scheme (space-time code) for the 2x2 MIMO cannel which 
allows for a very simple detector. Two complex symbols s\ and s 2 are transmitted in the first 
transmission (an odd transmission) and in the second transmission (the next even transmission) 
these symbols are more or less repeated. More precisely 
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The received signal is now 
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Rewriting this results in 
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(13) 



(14) 



1 Hcrc //t denotes the Hcrmitian transpose of H. It involves both transposition and complex conjugation. 
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or more compactly 

S\a + S26 + n, with 

{yu,y2i,yl 2 ,y 22 ) T , 

(hii,h2i,h* 2 , h* 22 ) 1 ', 
(hi2,h 2 2,-h* 11 ,-h* 21 ) T , and 

(nil, nai, 71*2, 7i2 2 ) r - ( 15 ) 

Since a and b are orthogonal the symbol estimates si and s" 2 can be determined by simply slicing 
(gfy)/(gJa) and (y^&)/(6 6) respectively. 

Another advantage of the Alamouti method is that the densities of g)a and 6^6 are (identical 
and) chi-square with 8 degrees of freedom. This results in a diversity order 4, i.e. 

Pr{(O a ) ^ (Si,5 2 )} « 7" ■ SNR -4 , (16) 

for fixed rate and large enough SNR. 

A disadvantage of the Alamouti method is that only two complex symbols are transmitted every 
two transmissions, but more-importantly that the symbols transmitted in the second transmission 
are more or less repetitions of the symbols in the first transmission. Section [T] however suggests 
that we can improve upon ordinary repetition. 
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4 The rotated and scaled Alamouti method 

4.1 Method description 

Having seen in section [1] that scaled-repetition improves upon ordinary repetition in the SISO 
case, we use this concept to improve upon the standard Alamouti scheme for MIMO transmission. 
Instead of just repeating the symbols in the second transmission we scale them. More precisely, 
when si and S2 are elements of -4i6-QAM = {a + jb\a £ A4-pAM,b G A^-pam}, we could transmit 
for some value of the signals 
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(17) 



where M 2 (ct) — 2a — D 2 (a) with D 2 (a) = 5(3 when j3 is the complex sign of a. 

A first question is to determine a good value for 8. Therefore we determine for < 9 < tt/2 
the minimum modulus of the determinant mindet(#) 

mindet(0) = min | det(X(si, s 2 , 6) - X(s[, s 2 , 6))\, (18) 

(si,s 2 ),(si,4) 

where X = ( Xl1 Xl2 ) is the code matrix. The minimum modulus of the determinant as a 



X21 X22 

function of can be found in Fig. O The maximum value of the minimum determinant (i.e. 7.613) 
occurs for 

opt . = 1.028. (19) 
We will use this value for in what follows. 
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Minimum Determinant, 16QAM-codesymbols, R=4 bits/transm. 




Figure 5: Minimum modulus of the determinant for rotated and scaled Alamouti as a function of 
8 horizontally. 



4.2 Hard-decision Performance 

We have compared the message-error-rate for several R = 4 space-time codes in Fig. [Sj By 
message-error-rate we mean the probability Pr{X ^ X}. Note that for each "test" we generate a 
new message (8-bit) and a new channel matrix. The decoder is optimal for all codes, it performs 
ML-decoding (exhaustive search). The methods that we have considered are: 

1. Uncoded, in green. We transmit 
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(20) 



where in, x±2, X21, and £22 are symbols from .Ai-qam- 

2. Alamouti, in blue, sec (| 1 2[> . where si and S2 are symbols from Ai6-qam- 

3. Tilted QAM, in cyan. Proposed by Yao and Wornell [7]. Let s a ,Sb,s c , and Sd symbols 
from »44-qam ■ Then we transmit 
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(21) 



for #1 = \ arctan(i) and 82 = h arctan(2). 



4. Rotated and scaled Alamouti, in red, see p7|) for 9 = 1.028, and with si and s 2 from 

^4l6-QAM- 



5. Golden code, in magenta. Proposed by Belfiore et al. [2]. Now 
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Figure 6: Message error rate for several R=4 space-time codes. 



6. Telatar, in black. This is the probability that the Telatar capacity of the channel is smaller 
than 4. 

Clearly it follows from Fig. [5] that the winner is the Golden code. However rotated and scaled 
Alamouti is only slightly worse, roughly 0.2 dB. Important is that Alamouti coding is roughly 2 
dB worse than the Golden code. 



5 Decoding complexity 

Clearly the Golden code is better than rotated and scaled Alamouti. However the Golden code 
in principle requires the decoder to check all 256 alternative codewords. Here we will investigate 
the complexity and performance of a suboptimal rotated and scaled Alamouti decoder. Denote 

6 = exp(jf3 opt .). 

A. In the rotated and scaled Alamouti case the received vector is 
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We can write this as 
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For the cos(0) of the angle between a and b we can write 
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(24) 



B. Instead of decoding (si, S2) we can also decode (ti,t 2 ) 
to (si, S2). Therefore we rewrite (fT7|) and obtain 



(M 2 (si), M 2 (s 2 )) which is equivalent 
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We can write this as 
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and for the cos(</>') of the angle between a' and &' we can write 
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C ° h{0 } ~ ^h.tf + A\h 21 \> + \h 12 \* + \h 22 \f [Z<) 
C. It now follows from the inequality 2r!r 2 < r\ + r 2 (where T\ and r 2 are reals), that 
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If 

\h 12 \ 2 + \h 22 \ 2 > |/in| 2 + |/i2i| 2 , (29) 

then cos(» < 21^zii = 0.393, else cos(<//) < = 0.393. Therefore it makes sense to decode 

(si,s 2 ) when (|2"9")) holds and (ti,t 2 ) when (f2"9")l does not hold. Using zero-forcing to decode, the 
noise enhancement is then at most 1/(1 — 0.393 2 ) = 1.183 which is 0.729 dB. We shall see later 
that noise enhancement turns out to be un-noticeable in practise. 

D. The decoding procedure is straightforward. Focus on the case where we decode (si, s 2 ) for 
a moment. For all 16 alternatives of (D 2 (si), D 2 (s 2 )) the vector 

z = y + D 2 (si)c + D 2 (s 2 )d = sia + s 2 b + n (30) 

and is determined. Then compute the sufficient statistic 
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Next use the inverted matrix M = ( ~ 7 ^ f ~ ^ /D where D = (a^a) (b}b) — (b^a) (g^b) to 
obtain [ S } J = M ( ) . Next both si and s 2 are sliced under the restriction that only 



s 2 J \ b'z 

alternatives that match the assumed values D 2 (s\) and D 2 (s 2 ) are possible outcomes. This is 
done for all 16 alternatives (D 2 (s\), D 2 (s 2 )). The best result in terms of Euclidean distance is 
now chosen. 

In considering all alternatives (D 2 (si), D 2 (s 2 )) we only need to slice when the length of z — 
s\a — s 2 b is smaller than the closest distance we have observed so far. This reduces the number 
of slicing steps. We call this approach METHOD 1. 

E. The number of slicing steps can even be further decreased if we start slicing with the 
most promising alternative (D 2 (s±), D 2 (s 2 )). This approach is called METHOD 2. Therefore we 

note that the "direct" s\ -signal-component in X is [ Sl * / )• Therefore we can slice 



-s\/2 

(e\y) / (e\§.i) m order to find a good guess for D 2 (si). Similarly we slice {ely)/(ele 2 ) to find a 
good first guess for D 2 (s 2 ). Here 

e, = (/nie,/i 2 ie,-r i2 /2,-/i* 2 /2) T , (32) 
e 2 = (-h 12 /2,-h 22 /2,-h* n ,-h* 21 ) T . (33) 

Then we consider the other 15 alternatives and only slice if necessary. Note that similar methods 
apply if we want to decode {t\,t 2 ). 

F. We have carried out simulations, first to find out what the degradation of the suboptimal 
decoders according to method 1 and method2 is relative to ML-decoding. The result is shown in 
Fig. E Conclusion is that the suboptimal decoders do not demonstrate a performance degradation. 
We have also considered the number of slicings for both method 1 and method 2. This is shown in 
Fig. [5) It can be observed that method 1 leads to roughly 7 slicings (as opposed to 16). Method 
1 further decreases the number of slicing to roughly 3.5. 



6 Conclusion 

Rotated and scaled Alamouti has a hard-decision performance which is only slightly worse than 
that of the Golden code, but can be decoded with an acceptable complexity. We finally remark 
that we have obtained similar results for codes based on mapping M 3 (-) for 9-PAM. 
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MER, Rot. Seal. Alam., 16QAM code-symbols, R=4 bits/transm., 1000 errors 
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Figure 7: Message error rate for three Rotated Scaled Alamouti decoders (R = 4), horizontally 
SNR. 



Av.nr. slicings, Rot. Seal. Alam., 16QAM code-symbols, R=4 bits/transm., 1000 errors 
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Figure 8: Number of slicings for two Rotated Sclaed Alamouti decoders (R = 4), horizontally 
SNR. 
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